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ABSTRACT 

Radial velocities of 40 intracluster planetary nebulae (ICPNe) in the Virgo cluster were obtained with the new 
multi-fiber FLAMES spectrograph on UT2 at VLT. For the first time, the A 4959 A line of the [OIII] doublet 
is seen in a large fraction (50%) of ICPNe spectra, and a large fraction of the photometric candidates with 
m(5007) < 27.2 is spectroscopically confirmed. 

ICPNe with the velocity dispersion of the Virgo cluster are found in our CORE field l.°2 from M87. These 
may have originated from tidal mass loss of smaller galaxies in the M87 subcluster halo. In a field 0.25 deg 
from M87, we see an extended stellar halo of M87 in approximate dynamical equilibrium, but with few ICPNe. 
Finally, in a field near M84/M86, the ICPNe velocities are highly correlated with the galaxy velocities, showing 
that any well-mixed intracluster population is yet to form. Overall, the measured velocity distributions confirm 
the non-uniform dynamical structure and on-going assembly of the Virgo cluster. 

Subject headings: (ISM:) planetary nebulae: general; galaxies: cluster: general; galaxies: cluster: individual 
(Virgo cluster); galaxies: evolution 



1. INTRODUCTION 

After the serendipitous discovery of intracluster planetary 
nebulae (ICPNe) in the Virgo (Arnaboldi et al. 1996) and For- 
nax (Theuns & Warren 1997) clusters, our group embarked on 
a systematic study of the diffuse stellar population in the Virgo 
cluster. The aim of this research is to measure the fraction of 
cluster stars between galaxies, their radial distribution, and 
their motions. The discovery of ICPNe, the subsequent nar- 
row band imaging surveys to find ICPN samples (Feldmeier et 
al. 1998, 2003a; Arnaboldi et al. 2002, 2003; Okamura et al. 
2002), and the HST observations of Virgo cluster empty fields 
(Ferguson et al. 1998; Durrell et al. 2002) have revitalized the 
study of diffuse intracluster light (ICL). The preliminary re- 
sults on the systematics of the ICL in different environments, 
from loose groups to rich clusters (Gonzalez et al. 2000, Gal- 
Yam et al. 2003, Castro-Rodriguez et al. 2003, Feldmeier et al. 
2003b), and the predictions from recent cosmological N-body 
and hydrodynamical simulations (Napolitano et al. 2003, Mu- 
rante et al. 2004, Sommer-Larsen et al. 2004, Willman et al. 
2004) have proven the study of the ICL as a valuable tool to 
investigate galaxy and galaxy cluster evolution. 

ICPNe are the only component of the ICL whose kinemat- 
ics can be measured at this time. This is important since the 
high-resolution N-body and hydrodynamical simulations pre- 
dict that the ICL is unrelaxed, showing significant substruc- 
ture in its spatial and velocity distributions in clusters similar 
to Virgo. While spatial structures have been observed in the 
ICPN number density distribution both in a single field (Oka- 
mura et al. 2002) and as field-to-field variations (Aguerri et al. 

* BASED ON DATA COLLECTED WITH THE FLAMES SPECTRO- 
GRAPH AT THE UT2 OF THE VLT AT CERRO PARANAL, CHILE, OP- 
ERATED BY ESO, DURING OBSERVING RUN 71.B-0147(A) 



2004), substructure in velocity space still needs to be investi- 
gated. 

Early attempts by Freeman et al. (2000) with the AAT and 
2dF spectrograph provided only a few spectroscopically con- 
firmed emission line objects in the Virgo cluster core. Their 
identification as ICPNe relies on the presence of the weaker 
[OIII] line (A 4959A) in the summed spectrum of all the de- 
tected single-line candidates. Because modern emission line 
candidate samples in Virgo contain a modest fraction of Lya 
emitters at z = 3.1 (Aguerri et al. 2004), spectral identification 
of ICPN either needs high spectral resolution to resolve the 
typical broad, asymmetric Lya lines, and/or detection of the 
[OIII] A 4959/5007 A doublet 1 . 

ICPN samples are sparse, with only a few tens of ICPNe 
in a 0.25 deg 2 field, and their fluxes are faint, from 1 x 10~ 16 
to 5 x 10" 17 erg cm" 2 s" 1 in the [OIII] 5007 A line. Spectro- 
scopic observations thus require 8 meter class telescopes and 
spectrographs with a large field-of-view (FOV), of a fraction 
of a square degree. The [OIII] emission lines from PNe are 
only a few km s" 1 wide, so relatively high resolution spectra 
(R = 7000 to 10,000) are desirable to reduce the sky contam- 
ination. The FLAMES -GIRAFFE spectrograph on VLT, with 
its medium-high spectral resolution, its FOV of 25 arcmin in 
diameter and 130 fibers in the MEDUSA mode, is therefore 
very well suited to this project. Here we present the results of 
the spectroscopic follow-up with FLAMES of the ICPN can- 
didates selected from three survey fields in the Virgo cluster 
core. 

1 The emission line sample of Kudritzki et al. (2000) was even dominated 
by Lya emitters; however, the luminosity function (LF) of this sample fol- 
lows the LF of Lya emitters in blank field surveys, not that of ICPNe in Virgo 
(Castro-Rodriguez et al. 2003). 
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Table 1 

Observed fields and spectroscopic confirmation rates 



FCJ a CORE" SUB" 



q(J2000) 12:30:47.6 12:27:47.4 12:25:33.8 

<5(J2000) +12:38:32.4 +13:21:40.1 +12:47:01.2 



m,, m (5007) b 


27.0 


27.2 


28.1 


W Fi b[< 27.2] b 


18 


34 


18 


N PN b 


15 


12 


13 


Af PN (A4959) b 


10 


5 


4 


N Plb (< m min f 


13 


34 


18 


Nm(< m min ) b 


11 


11 


13 


Rate(< mminf 


84% 


32% c 


72% d 


Predicted Rate(< m m ,„) b 


80% 


17% 


100% d 



"References for photometric catalogues: Arnaboldi et al. (2002, FCJ), 
Aguerri et al. (2004, CORE), Okamura et al. (2002, SUB). 

b Symbols denote: limiting magnitude of photometric sample, m/ im (5007); 
number of fibers allocated to objects with m < 27.2 and not in bad CCD re- 
gions, WFib[< 27.2]; number of detected PNe, Npn; number of detected PNe 
with double lines, ^Vpn(A4959); number of fibers allocated to objects brighter 
than min[27.2, m/, m (5007)], A , Fib(< fmin)', number of detected such PNe, 
A'pn(< m,„i„)', spectroscopic detection rate to m m ,„, Rate(< m m ,„) = /Vpn(< 
mmin)/N?\b(< Wmfn); corresponding rate predicted from photometry using 
simulations, Predicted Rate(< m m ,„) = number of PNe in statistically decon- 
taminated sample divided by total number of candidates in the field (Aguerri 
et al. 2004). 

c Higher than the predicted rate because in the fiber allocation priority was 
given to the brighter objects in the catalogue. 

d All candidates are true ICPNe, based on detections in both [OIH] and Ha 
narrow band images. Detected rate lower because of lower S/N. 



2. OBSERVATIONS 

Table 1 gives the three selected fields in the Virgo core, 
referred to as FCJ, CORE, and SUB (Aguerri et al. 2004). 
Spectra were acquired in service mode (10 hrs were al- 
located to this observing run, 71.B-0147, in priority A), 
with the FLAMES spectrograph at UT2 on VLT in the GI- 
RAFFE+MEDUSA configuration. We used the low resolu- 
tion grism LR 479.7, covering a wavelength range of 500 A, 
centred on 4797 A, and a spectral resolution of 7500. This 
gives a velocity resolution of 40 km s" 1 , and a typical velocity 
error of 12 km s . The redshifted [OIII] emissions of ICPNe 
in the Virgo cluster core fall near the red edge of the grism 
response. 

The FLAMES FOV covers the FCJ field entirely, and a sig- 
nificant fraction of both CORE and SUB. The total observ- 
ing time was 2.5 hrs for FCJ, and 2.6 hrs each for CORE 
and SUB; the exposure times were based on the S/N esti- 
mate for detecting the [OIII] 5007 A line flux of 4.2 x 10" 17 
erg cm" 2 s~\ i.e. m(5007) = 27.2, with a S/N~ 5. Given 
that the [OIII] 4959/5007 A emission lines have relative in- 
tensities 1:3, we expect to be able to detect the [OIII] 4959 
A emission for the brighter candidates only. The data reduc- 
tion was carried out with the GIRAFFE pipeline, for CCD 
pre-reduction, fiber identification, wavelength calibration, ge- 
ometric distortion corrections, co-addition and extraction of 
the final ID-spectra. 

Comparing with earlier measurements by Freeman et al. 
(2000), the RMS velocity difference for four ICPNe in com- 
mon is 24 km s" 1 , where their velocity error was 40 km s" 1 . 

3. SPECTROSCOPIC RESULTS 

A total of 40 ICPN candidates were detected in the 
FLAMES spectra. In Fig. ^ we show single ICPN spectra 
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FIG. 1.— FLAMES spectra for 14 ICPNe observed in the FCJ (two left 
columns) and CORE fields. The spectrum in the lower right corner is a Lyo 
object, which shows a very broad line profile. The m(5007) magnitudes are 
marked on the individual frames. 



with the [OIII] 4959/5007 A doublet, and the resolved spec- 
trum of a Lya galaxy with a broad asymmetric line. 

In Table 1 we give results for the three pointings individu- 
ally. In the FCJ/CORE/SUB fields, we had 18/34/18 fibers 
allocated to sources brighter or equal to 27.2 and in good 
regions of the CCD. In total, we detected 15/12/13 sharp 
line emitters, and 0/2/0 Lya emitters which show one re- 
solved asymmetric line. The remaining spectra did not show 
any spectral features in the wavelength range covered by 
FLAMES. The fraction of confirmed spectra with both com- 
ponents of the [OIII] doublet detected is 67%/41%/18%. 

The extracted ID spectra for the CORE and SUB fields are 
noisier than those for the FCJ, most probably because the se- 
lected reference stars for the FLAMES MEDUSA configura- 
tion were not in the astrometric system of the ICPN candi- 
dates, therefore fibers were not optimally positioned. For the 
SUB field, also the observing conditions were slightly worse 
than for the other fields, resulting in a low fraction of con- 
firmed spectra with detected [OIII] doublet. In the CORE 
field we checked that the summed spectrum of all sharp line 
emitters also shows the A 4959 A emission in addition to the 
A 5007 A line. The line ratio is 3.5, but with large error be- 
cause of the noise in the summed spectrum: these spectra are 
compatible with being all ICPNe. 

Table 1 also gives the spectroscopic confirmation rates in 
the three fields, for candidates down to the brighter of 27.2 
or m/, m (5007) from the photometry. This rate varies strongly 
from field to field despite similar m/,„,(5007). The most im- 
portant reason for this is different contamination of the sam- 
ples by faint continuum stars, erroneously classified as ICPNe 
because of a shallow off-band image. Aguerri et al. (2004) 
investigated the contamination of the photometric samples 
caused by faint continuum stars, [Oil] emitters at z = 0.347, 
and high-z emitters, using simulations and blank field surveys. 
Table 1 shows the expected spectroscopic confirmation rates 
based on these simulations. These are in close agreement with 
the actual confirmation rates, showing that the photometric 
samples are well understood. A high confirmation rate (small 
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FIG. 2. — Pointings surveyed in the Virgo core region. The field shown 
is 2° X 2° with M87 near FCJ and M86/M84 near SUB. The circles show 
the FLAMES FOVs. The square in the FCJ pointing marks the area of the 
imaging surveyed field. 



contamination) can be achieved when the off-band image is 
sufficiently deep. See Aguerri et al. (2004) for further details. 

4. ICPN LINE-OF-SIGHT VELOCITY DISTRIBUTIONS 

With these data, we can now for the first time determine 
radial velocity distributions of ICPNe and use these to in- 
vestigate the dynamical state of the Virgo cluster. Fig. [5] 
shows an image of the Virgo cluster core with the positions 
of our FLAMES pointings. The heliocentric radial velocity 
histograms obtained from the spectra in these fields are dis- 
played in Fig. [3] Clearly the histograms for the three pointings 
are very different. 

In the FCJ field, the velocity distribution of the PNe is not 
consistent with a single Gaussian with either the fitted veloc- 
ity dispersion ctfcj = 573 km s" 1 , or the more canonical Virgo 
g = 800 km s" 1 , based on a % 2 -test. Instead, it is dominated by 
a narrow peak, with v p = 1276 ±71 km s" 1 and a p = 247 ± 52 
km s -1 , which we identify with the halo of M87 below. In 
addition, there are 3 outliers, 2 at low velocity. All three are 
unusually bright, in the bright fall-off of the M87 PN luminos- 
ity function (PNLF), and they shift the FCJ luminosity func- 
tion to a brighter cutoff (Arnaboldi et al. 2002). The ICL sur- 
face brightness associated with these 3 outliers, e.g. the likely 
true ICPNe in the FCJ field, is /ig — 30.7 mag arcsec -2 if we 
use the conversion for M31 from Ciardullo et al. (1989), a\#. 
This is comparable with the surface brightness measurements 
of Ferguson et al. (1998) and Durrell et al. (2002) from intra- 
cluster red giants in fields further away from M87. 

The M87 peak of the FCJ velocity distribution contains 12 
velocities, which are described well by a Gaussian accord- 
ing to KS and x 2 tests. Their luminosity function is consis- 
tent with the PNLF in the inner 4' of M87 (Ciardullo et al. 
1998), with a KS probability pks = 0.46. The average veloc- 
ity agrees with that of M87, v sys = 1307 km s" 1 . The value of 
Up is consistent with the stellar velocity dispersion profile ex- 
trapolated outwards from ~ 150" in Figure 5 of Romanowsky 
& Kochanek (2001) and falls in the range spanned by their 
dynamical models for the M87 stars (the center of FCJ is 
15/0 ~ 65 kpc from M87, for an assumed M87 distance of 15 
Mpc). To infer the surface brightness corresponding to the 12 
PNe in the M87 peak, we again use ot\p for M31, because of 
the colour gradient observed in the outer parts of M87 (Goud- 
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FIG. 3. — ICPN radial velocity distributions in the three pointings (FCJ, 
CORE, and SUB). In the FCJ panel, the blue dashed line shows a Gaussian 
with v ml i = 1276 km s and cr rai j = 247 km s . In the CORE panel, the green 
dashed line shows a Gaussian with v m[ / = 1436 km s and cr m[ j = 538 km s , 
for VC galaxies dE and dSO within 2° of M87 (from Binggeli et al. 1987). 
In the SUB sub, the overplotted dashed histogram shows the radial velocities 
from a TNG spectroscopic follow-up (Arnaboldi et al. 2003). The dashed 
red line shows a Gaussian with v ral i = 1079 km s and a ra j = 286 km s for 
the M84 peak. The overplotted dashed-dotted lines show the SUB-FLAMES 
spectra including those of HII regions, which have radial velocities in the 
M84 and NGC 4388 velocity ranges. 



frooij et al. 1994). The resulting /i# = 29.2 mag arcsec -2 falls 
on the average M87 surface brightness profile extrapolated 
from ~ 400". fig and a p are thus also consistent with the 
extrapolated dynamical model of Romanowsky & Kochanek 
(2001) for the M87 stars, in the distribution of dark matter 
inferred by them, which is also similar to that determined by 
Matsushita et al. (2002). By contrast, the available data and 
dynamical models show an approximately flat dispersion pro- 
file for the globular clusters in M87, at ace — 350 km s -1 , 
corresponding to a shallower radial distribution. The main re- 
sult from our measurement of a p is that M87 has a stellar halo 
in approximate dynamical equilibrium out to at least 65 kpc. 

In the CORE field, the distribution of ICPN line-of-sight 
(LOS) velocities is clearly broader than in the FCJ field and 
consistent with a Gaussian (p KS = 0.9). It has v c = 1491 ±290 
km s" 1 and cj c = 1000 ±210 km s _1 ; the median is 1791 km 
s -1 . The CORE field is in a region of Virgo devoid of bright 
galaxies, but contains 7 dwarfs, and 3 low luminosity E/S near 
its S/W borders. None of the confirmed ICPNe lies within a 
circle of three times half the major axis diameter of any of 
these galaxies, and there are no correlations of their veloci- 
ties with the velocities of the nearest galaxies where these are 
known. Thus in this field there is a clear IC stellar component. 

The mean velocity of the ICPN in this field is similar to 
that of 25 Virgo dE and dSO within 2° of M87, < v dE .M87 >= 
1436 ± 108 km s" 1 (Binggeli et al. 1987), and with that of 
93 dE and dSO Virgo members, < v^ E vi rgo >= 1 139 ± 67 km 
s" 1 (Binggeli et al. 1993). However, the velocity dispersion 
of these galaxies is smaller, c7dE,M87 = 538 ± 77 km s" 1 and 

CdE.virgo = 649 ± 48 km s_1 ■ 
The inferred luminosity from the ICPNe in the CORE 

field is 1.8 x 1O 9 L B , (Aguerri et al. 2004). This is about 
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three times the luminosity of all dwarf galaxies in this field, 
5.3 x 10 8 Lb.o, but an order of magnitude less than the lumi- 
nosities of the three low-luminosity E/S galaxies near the field 
borders. Using the results of Nulsen & Bohringer (1995) and 
Matsushita et al. (2002), we estimate the mass of the M87 
subcluster inside 310 kpc (the projected distance D of the 
CORE field from M87) as 4.2 x IO'^Mq, and compute a tidal 
parameter T for all these galaxies as the ratio of the mean 
mass density within D to the mean density of the galaxy. We 
find T = 0.01-0.06, independent of galaxy luminosity. Since 
T ~ D~ 2 , any of these galaxies whose orbit now comes closer 
to M87 than ~ 60 kpc would be subject to severe tidal mass 
loss. Thus the ICPN population in the CORE field could be 
debris from the tidal disruption of small galaxies on nearby 
orbits in the M87 halo. The relatively small number of ICPNe 
in the FCJ field at D = 65 kpc could then mean that most of 
the tidally disrupted galaxies did not orbit as deep into M87. 

In the SUB field the velocity distribution from FLAMES 
spectra is again different from CORE and FCJ. The KS test 
gives low probabilities, pks = 0.13 and pks = 0.03, that the 
SUB histogram could be drawn from the velocity histogram 
in the CORE field or the Gaussian fitted to this, respectively. 
Instead, the SUB histogram of LOS velocities shows sub- 
structures that are highly correlated with the systemic ve- 
locities of M86, M84 and NGC 4388. The association with 
the three galaxies is strengthened when we plot the LOS ve- 
locities of 4 HII regions (see Gerhard et al. 2002) detected 
with FLAMES in this pointing. The highest peak in the dis- 
tribution coincides with M84, and even more so when we 
add the LOS velocities obtained previously at the TNG (Arn- 
aboldi et al. 2003). The 10 TNG velocities (2 velocities were 
measured in addition to the Arnaboldi et al. (2003) sam- 
ple) give vm84 = 1079 ±103 km s" 1 and cr M 84 = 325 ± 75 km 
s" 1 within a square of 4R e x 4R e of the M84 center. The 8 
FLAMES velocities in the M84 subpeak give v M 84 = 891 ± 74 
km s" 1 and <7m84 = 208 ± 54 km s" 1 , going out to larger radii. 
Note that this includes three over-luminous PNe not attributed 
to M84 previously. The combined sample of 18 velocities 
gives vm84 = 995 ± 69 km s" 1 and cr M 84 = 293 ± 50 km s" 1 . 
Most likely, all these PNe belong to a very extended envelope 
around M84 (see the deep image in Arnaboldi et al. 1996). 
It is possible that the somewhat low velocity with respect to 
M84 may be a sign of tidal stripping by M86, or of a recent 



merger with a smaller galaxy. In the latter case, the overlu- 
minous PNe might be due to a younger or a more metal-poor 
population (Dopita et al. 1992). 

5. CONCLUSIONS 

We have presented the first measurements of the veloc- 
ity distribution of ICPNe in three fields of the Virgo cluster. 
Overall, these velocity measurements confirm the view that 
Virgo is a highly non-uniform and unrelaxed galaxy cluster, 
consisting of several subunits that have not yet have had time 
to come to equilibrium in a common gravitational potential. 

A well-mixed IC stellar population is seen clearly only in 
the CORE field, in the outer parts of the M87 subcluster. 
Here the velocity distribution is consistent with a single clus- 
ter Gaussian, and the ICPNe might well have their origin in 
the tidal effects of the halo of this subcluster on its galaxy 
population. In the SUB field near M84 and M86, the ICPNe 
do not appear virialized; their velocities are highly correlated 
with those of the large galaxies in the field. In fact, there are 
regions in Virgo where no ICPNe are found (the LPC field of 
Aguerri et al. 2004). 

The measurements have also shown that M87 has a very 
extended envelope in approximate dynamical equilibrium, 
reaching out to at least 65 kpc. The cluster ICPN popula- 
tion in our FCJ field has density comparable with that in other 
fields further out, indicating a shallow ICL density distribu- 
tion. 



MA. would like to thank L. Pasquini and F. Primas for 
their help during the VLT-FLAMES observing block prepa- 
ration and service observing, and C. Cacciari, E. Pancino and 
E. Rossetti for their help with the GIRAFFE pipeline. We 
thank B. Binggeli and M. Capaccioli for useful discussions, B. 
Binggeli for providing his Virgo cluster galaxy catalogue with 
the dE & dSO radial velocity data, and an anonymous referee 
for constructive comments. We acknowledge financial sup- 
port by INAF - Projects of National Interests (P.I. MA), SNF 
grant 200020-101766, and Spanish DGES grant AYA2001- 
3939. NRN is supported by a MC fellowship within the EC 
Fifth Framework Program. This research has made use of the 
NASA/IPAC Extragalactic Database (NED). 



REFERENCES 



Aguerri, J.A.L., Gerhard, O.E., Arnaboldi, M., Napolitano, N.R., Castro- 
Rodriguez, N, & Freeman, K.C. 2004, AJ, submitted 
Arnaboldi, M. et al. 2003, AJ, 125, 514 
Arnaboldi, M. et al. 2002, AJ, 123, 760 
Arnaboldi, M. et al. 1996, ApJ, 472, 145 

Binggeli, B., Popescu, C. C, & Tammann, G. A., 1993, A&A, 98, 275 
Binggeli, B., Tammann, G.A., & Sandage, A. 1987, AJ, 94, 251 
Castro-Rodriguez, N., Aguerri, J.A.L., Arnaboldi, M., Gerhard, O.E., 

Freeman, K.C, Napolitano, N.R., & Capaccioli, M. 2003, A&A, 405, 803 
Ciardullo, R., Jacoby, G. H., Feldmeier, J. J., & Bartlett, R. E. 1998, ApJ, 

492, 62 

Ciardullo, R., Jacoby, G. H., & Ford, H. C. 1989, ApJ, 344, 715 

Durrell, P.R., Ciardullo, R. Feldmeier, J. J., Jacoby, G.H. & Sigurdsson, S. 

2002, ApJ, 570, 119 
Dopita, M. A., Jacoby, G. H., & Vassiliadis, E. 1992, ApJ, 389, 27 
Feldmeier, J.J., Ciardullo, R., Jacoby, G.H. & Durrell, P.R. 2003a, ApJS, 145, 

65 

Feldmeier, J.J., Ciardullo, R., Jacoby, G.H. & Durrell, P.R. 2003b, in IAU 

Symposium 217, Recycling Intergalactic and Interstellar Matter, in press 
Feldmeier, J.J., Ciardullo, R. & Jacoby, G.H. 1998, ApJ, 503, 109 
Ferguson, H.C., Tanvir, N.R., & von Hippel, T. 1998, Nature, 391, 461 
Freeman, K.C. et al. 2000, in ASP Conf. Ser. 197: Dynamics of Galaxies: 
from Early Universe to the Present, 389 



Gal- Yam, A.H. et al. 2000, ApJ, 536, 561 

Gerhard, O., Arnaboldi, M., Freeman, K.C, Okamura, S. 2002, ApJ, 580, 
L121 

Gonzalez, A. H, Zabludoff, A. I., Zaritsky, D., & Dalcanton, J. J. 2000, ApJ, 
536, 561 

Goudfrooij, P., Hansen, L., Jorgensen, H.E. & Norgaard-Nielsen, H.U. 1994, 

A&AS, 105, 341 
Kudritzki, R.-P, et al. 2000, ApJ, 536, 19 

Matsushita, K., Belsole, E., Finoguenov, A., Bohringer, H, 2002, A&A, 386, 
77 

Murante, G. et al. 2004, ApJ, 607, L83 
Napolitano, N R. et al. 2003, ApJ, 594, 172 
Nulsen, P. E. J., & Bohringer, H, 1995, MNRAS, 274, 1093 
Okamura, S. et al. 2002, PASJ, 54, 883 
Romanowsky, A. J., & Kochanek, C. S. 2001, ApJ, 553, 722 
Sommer-Larsen, J., Romeo, A.D. & Portinari, L. 2004, MNRAS, submitted 
i astro-ph/0403282 i 

Theuns, T. & Warren, S.J. 1997, MNRAS, 28 4, LI 1 

Willman, B. et al. 2004, MNRAS, submitted i astro-ph/0405094 i 



